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As shown in previous work, number projection can be carried out analytically for states defined 
in a quasi-particle scheme when the states are expressed in a coherent state representation. The 
wave functions of number-projected states are well-known in the theory of orthogonal polynomials 
as Schur functions. Moreover, the functions needed in pairing theory are a particularly simple class 
of Schur functions that are easily constructed by means of recursion relations. It is shown that 
complete sets of states can be projected from corresponding quasi-particle states and that such 
states retain many of the properties of the quasi-particle states from which they derive. It is also 
shown that number projection can be used to construct a complete set of orthogonal states classified 
by generalized seniority for any nucleus. 



I. INTRODUCTION 



The loss of particle number as a good quantum number is a serious deficiency of the BCS and HFB (Hartree-Fock- 
Bogolyubov) approximations in their applications to finite nuclei. However, it can be restored with relative ease in a 
coherent-state representation \V^ in a manner that preserves much of the elegance and simplicity of the models. 
Following early suggestions by Mottelson 0], number conserving extensions of BCS and HFB theory have been 
^^ ' considered by many authors: cf., for example, the number-projected BCS approximation of Dietrich, Mang, and 
^D . Pradal ^, the projected quasi-particle model of Lande, Ottaviani and Savoia g, the antisymmetrized geminal power 
CTN ' states of Coleman et al. 0, the broken pair model of Lorazo, Gambhir, and others RB], Talmi's generalized seniority 
scheme |9|, and the coherent correlated pair method of Vary and Plastino ||l^. A review of the subject has been given 
by AUaart et al. ^. 

A common feature of these schemes is to approximate the ground state, for a system in which pair correlations 
^ • dominate, by a state of the form 



|A/-)=H^|a:), (1) 



where \x) is a quasi-particle vacuum state and H^r is a projector onto the space of states of particle number N . We 
show herein that such a projection can be applied analytically, when the wave functions are expressed in a coherent- 
state representation. More importantly, projection from quasi-particle states gives a complete set of iV-particle states 
which reflect many of the properties of the states from which they were projected. Thus, a potentially exact formalism 
emerges which is able to take advantage of the physical insights gained from the BCS and HFB approximations. 
d ' Analytic number projection when combined with analytic methods for projecting out states of good angular momen- 

tum from deformed intrinsic states jl^ raises the possibility of developing good theories of nuclei which exhibit both 
rotational and superconducting properties. It also suggests a way of extending the physically significant definition of 
generalized seniority to the classification of complete sets of orthonormal states. 

For simplicity, we restrict consideration in this paper to the BCS model. However, the methods extend to the 
general HB theory by the number projection methods given in ref. ||lj]. 

II. QUASISPIN ALGEBRAS AND THEIR COHERENT STATE REPRESENTATION 

For each j-shell, there is an SU(2) quasispin algebra spanned by the operators 

'^+ — 2^m=l/2"im"im' '^- — Z^m=l/2 " " ' 

^0 = ^ E («]»«"" - °'"«]™) ' (2) 

m=-j 

where 



are creation and annihilation operators for a nucleon (e.g., neutron) which satisfy the fermion anticommutation 
relations 

{a-'",a]j = <5„„. (4) 

Note that we follow Einstein's convention of using upper indices to label the components of a tensor eontragredient 
to a tensor labeled by lower indices. With this convention, the annihilation operator a^™ with superscripts is the 
Hermitian adjoint of the creation operator aj„ and the annihilation operator ajm with subscripts is the m component 
of a spherical tensor. The contraction of a pair of upper and lower indices is then a scalar; e.g., 

is the number operator for the j-sheU. The quasispin operators satisfy the usual SU(2) commutation relations 

[S'_^,Sl]^2Si, [SlSi]^±Si, (6) 

and the nucleon operators transform under SU(2) as components of quasispin-1/2 tensors: 

[o_, 0,jm\ — '^jm, [-3+ , (ljm\ = 1 



i„ t 1 _ l„t rcJ' . 1 _ 1 



[oQ,a „jj — 2'^im ' [^'1 



'jm\ ~ 2 jni ' L*-'0'"'J™J ~ 2 J™ 



a," 



(7) 



Representations of a quasispin algebra are characterized by highest and lowest weights. For example, the zero- 
particle state |0) is a lowest weight state for all the quasispin algebras; 

Si\0) = 0, ^^|0) = -s,|0), V.7, (8) 

where Sj = (2j + l)/4. Such a state is said to have multishell seniority (0, 0, . . .). More generally, a state {I'JM) of i' 
particles that satisfies the equations 

Si\iyJM)^0, Sl\vJM)^{-Sj+Vj/2)\vJM), Vj , (9) 

is said to have seniority v and multishell seniority [vj-^ jVj^,...), where v = '^jVj- The multishell seniority defines the 
values of the lowest weights and provides useful identifiers of the corresponding quasispin irreps. 
A coherent state for the combined quasispin algebras is defined by 

\xvpJM)^e^+^'''^\vpJM) with S+{x) ^^XjSl^ , (10) 

j 

where \vpJM) is a lowest weight state and p provides whatever labels are needed, in addition to vJM , to characterize 
the state. A coherent state is a quasispin transform of a lowest weight state. The simplicity of the BCS and HFB 
theories stems from the fact that the quasi-particle operators of the Bogolyubov-Valatin transformation are also 
quasispin transforms of nucleon creation and annihilation operators. In particular, a coherent state of the seniority- 
zero irrep 

|.T)=e^+(")|0) (11) 

is the vacuum state of the quasi-particles defined, for real values of {xj}^ by the Bogolyubov-Valatin transformation 

Oijm — Ujajm + 'L'jO'jm ■ 

with V?, -\- v^ — 1 and Vj ~ XjUj. This well-known result follows from the observation that 

aj™e^+(^) = e^+(^)("jajm + Vja]^ - XjU^a]^) = u^e'^+^^'aj™ . (13) 

Although the transformation e^+^^' is not unitary, it is simply related to the unitary transformation 



g(a;) = e^+(«-^-(«, (14) 

for which 

1 a; 

cosPj = —j^=^= = Uj , sin/3j — '' = Vj . (15) 

The relationship is given by the expansion 

g(x) = e^+(") \[{l + x]fo e-^-(-) . (16) 

3 

Thus, the transformation 

|0) -. \x) - e^+(-)|0) = g(x)|0) n(l + ^^Y' (17) 

j 

is unitary, to within a factor ■\/$^(a?^, where 

$^(x2) = (xia;) = [|(1 + x|)2«^- (18) 

and x^ ~ (a;^^,^^^, . . .). Likewise, the BV transformation is seen as the unitary SU(2) transformation 



^jm ^^ '^jm ~ 



Kx)a]^g^{x) , aj„i -^ ajm = g{x)aj„,g\x) . (19) 

In a VCS (vector coherent state) representation |l^jlj], a fowest weight state \vpJM) is regarded as an intrinsic 
state for an irrep of the combined quasi-spin algebras and represented by an intrinsic wave function ^^pjM- An 
arbitrary state \tp) is then represented by a hofomorphic wave function ^ defined over a set of complex coordinates 

z = {zj^ , Zj2 , . . .) with values 

*(^)= E C.pJM(i^pJM|e^-W|^). (20) 

upJM 

For example, if the zero-particle vacuum |0) has intrinsic wave function ^Qj an arbitrary quasi-particle vacuum state 
a;) — e'^+^^'jO) has wave function ^^ with values 

*x(2) = eo(0|e^-(^)e^+(")|0) = io-^^zx) , (21) 

where $*(za;) is the overlap 

^^{zx)^{z\x)=\{{l + z,x,f-^^ (22) 

j 

with zx = {zj-^Xj-^ , Zj^Xj^ ,•■•)• It follows that the norm of a quasi-particle vacuum state \x) is the value of the function 
^^ at x'^. 

In the VCS representation, an element X of an SU(2) quasispin algebra is represented as a differential operator 
r{X) defined by 



[r(X)*](z) = J2 ^.pJM{i^pJM\e^-^'^X\^) 

upJM 

= J2 ^•^pJm{i^pJM\{X + [S-{z),X] 



vpJM 



-i[5_(z),[5_(z),X]])e^-W|V). (23) 

Evaluation of the right hand side of this expression by means of the identities 



r{si,) = -s, + z, — , nsL)^ — , r{si)^z,[2s,-z~], (25) 



(i//9jM|S'^e^-(^)|V') =0, {lypJM\S^f^e^-^'''>\^P) = (-s^- + z/j/2)*(z) , 
gives the representation 

where Sj is a diagonal operator on intrinsic wave functions with eigenvalues 

Sj^iypJM = {-Sj + l^j/'2)S,i,p.JM ■ (26) 

The number operator fij is represented 

d_ 



T{nj)^0,+2z,— . (27) 



III. NUMBER-PROJECTED STATES 

Let \n{x)) denote the n-pair state 

|n(x))- [S+{x)]''\0)=n\n2n\x). (28) 

States of good nucleon number can also be projected from excited quasi-particle states. Consider first the one quasi- 



particle state al^la;). To within a factor Uj = 1/w 1 + x^, the (2n + l)-particle component of such a state is given 

by 

n2„+ia]„|.T) - n2„+ie^+(-)[z.2 + x,u,v,]a]jO) - 1 [5+(a;)]"|jm) , (29) 



where Ijm) = aj„|0). This is a special case of a general result. 

Claim 1: Let ZngjMi'^'') denote some combination of /i nucleon creation operators 

zlpjMi^^) - E Cjm" [a], ® a], ® • • • ® aljjM (30) 

j 

which creates a state 

|M/3JM)-Zt^^^(at)|0) (31) 

having the property that 

S^{uv)\^(3JM) ^0; (32) 

the subscript (3 provides whatever labels are needed in addition to jiJM to specify the state. Let Zj^sjMi'^'^ ) denote 
the corresponding quasi-particle operator in which each a- is replaced by Uja,^. Then 

^2n+pZl^jMi^a^)\x) =^ -,[S+(x)^\^lPJM) . (33) 



Proof: The claim follows from the identity 
which implies that 



t ,s+(.) ^ e^+(.)( t^ _ ^^.^^.^^.^) ^ ^S+(.)^-S_(«.)^t ,S_(™) (34) 



Q.E.D. 

It will be noted that to obtain complete sets of linearly independent states by number projection in this way, the 
states {\fi/3JM)} should include not only states that satisfy 

Si\fi(3JM) = yj , (36) 

but also states such as S+{y)\iJ, — 2,/3'JM} for y ^ x. Imposing the constraint 

S-{uv)S+{y)\[i-2,ii'3M)^Q (37) 

for all such states is one way, but not the only way, of ensuring that no linear combination of such states is of 
the form Sj^{x)\[i — 2,(3' JM). The latter states are not needed because the (2n + /i)-particle states projected from 
e'S'+(^)5'_i_(a;)|^ — 2,(3' JM) are proportional to those projected from e^+'^^^\^i — 2,(3' JM). 
The above results lead to the following claim. 

Claim 2: Let {Z]^njj^j{a^)\Q)} denote a complete set of orthonormal states, as defined in Claim 1, that satisfy the 
constraint S^{uv)Zujj^j{a^)\{)) — 0. Then the A'^-particle states number-projected from the corresponding quasi- 
particle states {^IajmC'-'^^)!^)} forni a (nonorthonormal) basis for the A^-particle nucleus. 

IV. VCS WAVE FUNCTIONS OF NUMBER-PROJECTED STATES 
The number-projected state \n{x)) has VCS wave function 

*n(.)(2) = ^o(0|e^-(^V(a^)) = £.oK{zx) , (38) 



where $fj = 7t,!P„$* and Pn^^ is the component of $" of degree n in its argument. It turns out, as shown in refs. [yg] 
that the function $^ is well-known in the theory of symmetric polynomials flSt] as a Schur function (or S-function) ; 
such functions are the characters of fully antisymmetric irreps of the unitary groups and easy to derive. 
The definition 



implies that 



where 



Kizx) ^ {Q\e^-^^\S+{x)r\Q) (39) 

K{^x) = [T{S+{x)Wn_,]{zx). (40) 



r(5+(x)) = Y. ^^nSi)) = E ^.^. (2-'^. - ^.^) • (41) 

J 3 ■' 



Thus, the needed Schur functions satisfy 

d_ 



j ■' 

Following the methods of ref. HB], this recursion relation is solved by making a change of variables from {zj} to 
the symmetric power functions 



E2s,z™. (43) 



3 

The recursion relation then becomes 



n-1 



$'* =: 



6^-^mC+iV^ 



$ 



71 — 1 ' 



where 



and has solutions 



y 



(44) 



(45) 



1>o = 1 



K? - ^2 



^2 '— KVlJ - V2 

$1 = (0f )3 - 30f,^^ + 2</>| , etc. 



(46) 



The general expression is obtained by means of the identity 

1 fi\ 
m [n ~ my. 

derived in ref. |l|. With this identity, the recursion relation is given in the useful form 

-i(^-l)! 



^^ = E(-i) 



[n — m) 



I Vm n—m 



(47) 



(48) 



and has explicit solution given by 



$f, = det 



1 



••• 



K (t>n-l 0n-2 0n-3 



Some properties of these functions are discussed in refs. ||l|,p|,[r5[ . 

It follows from Claim 1 that the VCS wave function for a one quasi-particle state Uja-Jx) has values 



From the observation that 



it follows that 



where 



This function has values 



where s' is the set of quasi-spins with 









$-'(z) = -^$^(z)=.$«(z), 

i + Zj 

s'k = Sk- ^Skj ■ 



(49) 

;iven by 

(50) 

(51) 
(52) 

(53) 

(54) 

(55) 
(56) 



6 



Thus, by considering $* to be one of a set of functions parameterized by values {sj = (2j + l)/4} of the quasispins, 
^^ is seen to be the member of the set with Sj replaced by Sj — 1/2. This substitution will be denoted by an operator 
A^ i.e., 

A^" : s, ^ s, - i , (57) 

so that 

$J = AJ$" . (58) 

The adjustment of the coherent state wave function <&* — > $■' by the replacement Sj — > Sj — 1/2 to take account of 
the occupation of one state of a single-particle level is an expression of the well-known blocking effect. 
The number-projected one quasi-particle state 

\n{x)ljm) = n!u,n2„+ia],Jx) = [S+{x)r a] JO) (59) 

is now seen to have VCS wave function with 

'^n(x)i3m{z) = n!^j„P„$J (za;) = ^j,n<^i{zx) , (60) 

where $^ is the Schur function 

K = ^'K ■ (61) 

Similarly, VCS wave functions are derived for all other number-projected quasi-particle states. For example, the 
seniority two states projected from a two quasi-particle state nluj^Uj^ [a]: (giaj ]jm\x) with J 7^ have wave functions 

*«(x)0u-2)./m(^) - [Oi «• ^niiM'^i'"' (zx) , (62) 

where 

$ili2 ^ ^il ^i2 $^ ^ $^" (g3) 

is the Schur function for the quasi-spin set s" with components 

4' = sfc-i(4ji+4jj. (64) 

The seniority-zero state S+{y)\n{x)) has VCS wave function given by 

^oY.y,z,{2s, - z, A)<i,^(^^) ^ ^^[^K^y) -Y^mMyin^rnWnii) , (65) 

where x = xz and y — yz. Thus, with the identity (p7|), the VCS wave function for the state S+{y)\n{x)) is given by 
UO\e'-^'^S+{y)\nix)) = ^o ^ (-1)"' , ""' ,, Myi"')'^n-m{i) ■ (66) 

m=0 ^ ' 

As expected, this wave function becomes identical to ^o^n+i when y = x. 

V. EVALUATION OF ENERGIES AND MATRIX ELEMENTS 

To illustrate the techniques, we consider the simple BCS Hamiltonian 

H = Y,£3nj-GS+S^, (67) 

j 

with 



E 






>5±-E^: 



(68) 



For given real values oi x — {xj-^ , Xj^ ,-■■), the quasi-particle operators are defined by 



ajm — Ujajm + 'VjO.jm i 



such that u^ + v'^ ~ 1 and Vj — XjUj so that 






l + x 



2 ' ^i'^i 



i + xy ^ 



i + xi 



(69) 



(70) 



A. The number-projected vacuum energy 

From standard BCS theory, the expectation of H for a quasi-particle vacuum state is given by 
{x\H\x) = [^2s,(2e, -G^;|)^;2 -G(^2s,«,«,)']<5^(a;2). 



With the above expressions for Uj and Vj , this expression becomes 



iSjXj \ 



$"(.t2) 



^^Asjejx']¥{x'^) - G^2sja;|$J'^(a 

3 3 

-G V'4siSjXiXj$'^(x^) . 

ij 

The energy of the n-pair state \n(x)), 



E^ix) 



{x\H\n{x)) 
{x\n{x)) 



(71) 



(72) 



(73) 



is now obtained by picking out the components of {x\H\x) and {x\x) of degree n in a;^. Recalling that P2n$'' — ^ti/^^-' 
we immediately obtain 



E-ix) 



$^(x2) 



^ [4s,s,x'^<^i^,{x^) G{n l)2s,x^<^^l,{x') 



-G^4s,s,x,x,$tf'_i(x2)j .) 



(74) 



Thus, the values of the x coefficients can be fixed such that the energy Eq — Eq{x) is minimized and the variational 
equation 



^(z|(i7-i?o")Kx)) 



0, Vj 



(75) 



is satisfied. 



B. Number-projected one quasi-particle energies 

As observed above, the (2?! + l)-particle component of the one quasi-particle state nlujaj^lx) is the state 
\n{x)ljm) — a'j^\n{x)). Thus, we consider the energies of number-projected one quasi-particle states defined by 

E" {n{x)ljm\n{x)ljm) ~ {n(x)ljm\(H ~ EQ)\n{x)ljm) (76) 

or equivalently by 

E- <^i{x') = {xW^{H - E-)a]jn{x)) . (77) 

Now the coherent state representation of the number operator fij ~ X^m '^jm'^"''"' when acting on a seniority zero 
state, is given by 

r(^,)-^2z,A. (78) 

Therefore 

{xW"'a],^{H-E-)\n{x)) = (l-^±.){z\{H-E-)\n{x)) 

\ ZSj OZj / z—x 

= (79) 

when X is assigned the value for which eqn. (|7^) is satisfied, ft follows that 

Ei <i>i{x') = {xy"^[H,alj\n{x)) . (80) 

If we now (temporarily) regard a; as a variable parameter, we obtain from Claim 1 the identity 

{xy"^[H,a]j\nix)) = u,{x\a^^-[H,a]j\n(x)) 

^nlP,,Uj{x\a^'^[H,alj\x), (81) 

where P„ picks out the component of degree n in x^ in what follows it. Note that for a particular value of x, the 
matrix element {x\a^"^[H,a,jj^]\x) is a number. However, by regarding x first as a variable and evaluating the matrix 
element as a function of x, it is meaningful to extract the component of this function of degree n; the rhs of eqn. (|8l| ) 
is then the value of this component at the assigned value of x. 

The quasi-particle vacuum properties of the state \x), which hold for arbitrary x, can now be used to put the matrix 
element on the right into the standard equations-of-motion form p6| 

{x\a'^[H,a]j\x) =. (x|{a-'"\ [H,a]j}\x) . (82) 

This matrix element is then evaluated by BCS methods to give 

uj (a;|{aJ™, [H, a]j}\x) = [e^w^ _ q^2^2 _^ q^^^^ ^ 2s,u,;w,] ^%x^) . (83) 

i 

The final step of number projection is now easy and gives the expression for the number-projected one quasi-particle 
energies 

E^<^i{x') = e,<^Ux') - nGx]^^^_,{x^) + nGY^2s,x,x,<^l_,{x') . (84) 

i 

This is an explicit and simple expression for _E" which involves only the values of known (Schur) functions at the 
specific X for which Eq{x) is a minimum. 



C. Matrix elements between seniority-two states 

Let {AJ,j j^^(at)|0}} be an orthonormal set of two-particle seniority-two states. The matrix elements between the 
corresponding number-projected quasi-particles states are given by 

= n\P,,(s,,,uuluf{x\H\x) + (xU^^'^'Hua), [H,Alj,,{uar)]]\x)) . (85) 

Thus, the matrix elements for seniority-two states are given by number projection of corresponding quasi-particle 
Tamm-Dancoff expressions, as defined in the equations-of- motion formalism pM. Moreover, the projections are 
accomplished by expanding the unprojected expressions in terms of Schur functions. 

D. Matrix elements between seniority-zero states 

After evaluating the overlaps {x\S-{x)S+{y)\n{x)) and {x\S-{y)S+{y')\n{x)), one can construct an orthonormal 
set of seniority-zero states for a 2n-particle nucleus. However, so that use can be made of Claim 1, it is convenient to 
start with operators {5-(-(j/)} that satisfy the condition 

S^{uv)S+{y)\0) = 0. (86) 

The overlap {x\S-{x)S+{y)\n{x)) = {n + l){x\iS^{y)\n{x)) is obtained from the value ^o(2;|5'+(j/)|?T-(a;)) at z = x of 
the coherent state wave function for the state S+{y)\n{x)); an expression for this wave function is given by eqn. (p9). 
The second overlap is given by 

{x\S^{y)S+{y')\nix)) = J] Y. %2/;'(^l«^""'«^"4™'4r..>(^)) 

j,m>0 j' ,7n'>0 

= J22sjyjy'^K^{x'). (87) 

3 

The matrix element {x\S-{x)HS^{y)\n{x)) = {n + l){x\H S^{y)\n{x)) is evaluated starting from the observation 
that 

{x\{H - ES)S+iy)\n{x)) = (n(a:)|(iJ - E^,)S+iy)\x) 

= {x\S^{y){H - E°^)\n{x)r 

and that 

d_ 



z\S-iy)iH - E-)\n{x)) ^Y.y^'^i'^'^ - '^-^)('\(H ^ ^o)Hx)) (89) 



vanishes when z = x because of eqn. (^3) . It follows that 

{x\S.ix)HS+iy)\n{x)} = (n + l)E^,^{x\S+iy)\n{x)} . (90) 

Writing 

{x\S.{y)HS+iy')\nix)) = n\P^{x\S-{y)HS+{y')\x) , (91) 

and defining 

S^iu^)^ J2 zi,'2/,a^"a^" , S+iuY) ^ E "'y>k«k , (92) 

j,rn>0 j,rn>0 

we can use Claim 1 to obtain 

{x\S^{y)HS+{y')\n{x)) = nlP,,{x\S-{u^y)HS+{u^y')\x) . (93) 

10 



Then, with the identity 

{x\~SAu^y)S+{u^y')H\x) ^Y.'^s,u^^y,v'j{x\H\x) , (94) 

we obtain 

{x\S-(y)HS+{y')\n[x)) = n\P,,[Y,2sjU^jyjy'j{x\H\x) 

3 

+ {A[S-{u^y)AHS+{uy)]]\x)\. (95) 

VI. GENERALIZED SENIORITY 

Talmi Q has defined the sequence of J = states proportional to {[5'+(x)]"|0); n = 1, 2, . . .} as states of generahzed 
seniority zero. Likewise, the sequences of J 7^ states {[S'+(a;)]"|2JM); n = 1,2,...}, where \2JM) is a two- 
particle state that is annihilated by all the quasispin lowering operators (i.e., S^_\2JM) — 0, Mj), are defined to be 
states of generalized seniority two. States of higher generalized seniority are similarly defined. Thus, the concept 
of generalized seniority identifies subsets of states of the corresponding seniority (more precisely summed multishell 
seniority v = '^aVj)- The generalized seniority zero states, for example, exclude states that are generated by 

combinations of the {5^_} operators that are not simply multiples of S+{x). Such states, sometimes described as 
broken S+{x) pair states M, are also described as generalized seniority-two or higher states, depending on how many 
S+{x) pairs are broken. The problem is that, while one can define two-particle J — Q parent states {\N — 2,J — 0)} 
as states that are annihilated by the S^{x) lowering operators, in parallel with lowest-weight quasispin states, the 
sequences of states {[S'+(a;)]"|2, J = 0);n = 1,2,...} are not orthogonal to the generahzed seniority-zero states. 
Consequently, as Talmi has emphasized, generalized seniority does not define a complete orthogonal scheme. 

Nevertheless, generalized seniority can be given a definition that does characterize a complete set of orthogonal 
states for each nucleus separately. For example, one can define a generalized seniority- two creation operator S'+(2/) 
for the 2ri-particle nucleus such that 

{x\S+{y)\n{x)) ^ Q . (96) 

With the identity 

{z\S+{p)[S+{x)T\Q) = n\PMS+{yp)\x) = nlPr,J2vj-^'^%^^) > (97) 

j ■' 

this equation reduces to the easy to solve equation for the yj coefficients 

Y,'2sjy,x,<^i{x')=0. (98) 

i 

States of higher generalized seniority may be defined similarly. 

Such a definition, has not been used to our knowledge probably because it would appear to destroy the simplicity 
and elegance of the concept. However, the facility to carry out number-projection analytically suggests that this may 
no longer be a significant concern. 

VII. APPLICATION TO A TWO-LEVEL MODEL 

To check the above methods, they have been applied to a simple two- level model having Hamiltonian H — "Y^- Sjhj — 

GS-i^S^ with £1 = and £2 = 1; the excitation energy £2 — £1 = 1 then sets the unit of the energy scale. For simplicity, 
we also set si = S2 = s and considered a model nucleus with n = 2s nucleons; this corresponds to a situation in 
which, in the G = limit, the lower level is fully occupied and the upper level empty. 
For such a two-level model, we can set 
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Xl = 1 , X2= X , 



(99) 



because replacing these parameters values by the substitution xi -^ k, X2 -^ kx would only change the overall norm 
of the state 



\n{x)) = [Sl + a;5'^]"|0) -> fc"[S'| + a;5^]"|0) 



(100) 



The energy Eq{x) of the state, with n — 2s, is then evaluated directly from eqn. ( |74| ) and the value of x, for which it 
is a minimum, determined. 

The structure of the state \n{x)) is revealed by expanding it on a basis of states labeled by the two-level quasispin 
quantum numbers 



{|to) — \s,s — m]s,m — s) , m 



,+s}. 



(101) 



The state \m) is the state with m nucleon pairs in the upper level and s — m pairs (or m hole pairs) in the lower level. 
The expansion 



i(a;)) = ^ Crn\m) 



(102) 



in this basis is easily inferred from the identity 

^^ix') ^ {xHx)} , 



(103) 



which implies that the squares of the Cm coefficients are given by the corresponding expansion of $|s(x^) as a 
polynomial in cc^™; 



'^2s{-')-^T.n 



Px^'". 



(104) 



The Cm coefficients are shown, for s = 7/2 and three values of G, in comparison with exactly computed wave functions 
for the model in fig. 111. 




FIG. 1. Ground state wave functions for n — 2s = lA particles in a two-level pairing model for three values of the pairing 
interaction. The expansion coefficients, as defined in the text, for normalized number-projected states are shown as small 
circles. The corresponding coefficients for exactly computed ground states are connected by continuous lines. 

For s — 7/2, the model corresponds to 14 nucleons in two j = 13/2 levels. Results computed for other values of 
n = 2s show that the number-projected wave functions rapidly become exact as s — » oo. 

As an indication of the improvement given by number projection over unprojected BCS results is given by comparing 
the ratio of the number of nucleons occupying the upper level to the number in the lower level. The mean number in 
the upper level is given for the above wave functions by 
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(n2)=^2m|a„|'a;''" 

^"|C„.P.^^ (105) 

and the mean number for the lower level is 

(m) =2s- (na) . (106) 

Fig. ?? shows the ratio {n2)/{ni) for several values of the pairing interaction in comparison to the exact ratios and 
those of the BCS approximation. As n = 2.s -^ cxd, the occupancies given by the BCS approximation are found 



G/Gcrit 

FIG. 2. The ratio of the mean number of particles in the upper level to the mean number in the lower level for the n — 2s — 14 
two-level pairing model defined in the text. The ratios for number-projected states are compared with exactly computed results 
and with those of the BCS approximation. 

to approach those of an exact calculation. However, for a finite number of particles, a considerable improvement is 
gained by number projection (as expected from previous results). 

VIII. CONCLUDING REMARKS 

The calculations that have been done, cf. for example refs. |g,|ll|, show clearly that number projection effects very 
substantial improvement to the accuracy of the BCS approximation in applications to finite nuclei. Moreover, the 
coherent state techniques introduced in refs. ||l| and [g were found to facilitate the practical steps of carrying out 
the number projection considerably. These techniques have been developed in this paper in the hope that they will 
be taken up by others in calculations of the low-energy spectra of nuclei in which pair-correlations are dominant. 
Of particular interest are the low-energy spectra of single-shell nuclei with both pairing and quadrupole interactions 
where the objective will be explore the emergence of deformed rotational states. 

Other recent developments that may be usefully deployed in concert with number projection are: the identification 
of a range of seniority-conserving interactions p7[ , and analytic techniques to carry out angular momentum projection 
from a class of deformed intrinsic states f|12| . 
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